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ABSTRACT 

We present an analysis of the gas physics at the base of jets from five T 
Tauri stars based on high angular resolution optical spectra, using the Hubble 
Space Telescope Imaging Spectrograph (HST/STIS). The spectra refer to a region 
within 100 AU of the star, i.e. where the collimation of the jet has just taken 
place. We form PV images of the line ratios to get a global picture of the 
flow excitation. We then apply a specialised diagnostic technique to find the 
electron density, ionisation fraction, electron temperature and total density. Our 
results are in the form of PV maps of the obtained quantities, in which the gas 
behaviour is resolved as a function of both radial velocity and distance from 
the jet axis. They highlight a number of interesting physical features of the 
jet collimation region, including regions of extremely high density, asymmetries 
with respect to the axis, and possible shock signatures. Finally, we estimate the 
jet mass and angular momentum outflow rates, both of which are fundamental 
parameters in constraining models of accretion/ejection structures, particularily 
if the parameters can be determined close to the jet footpoint. Comparing mass 
flow rates for cases where the latter is available in the literature (i.e. DG Tau, 
RW Aur and CW Tau) reveals a mass ejection-to-accretion ratio of 0.01 - 0.07. 
Finally, where possible (i.e. DG Tau and CW Tau), both mass and angular 
momentum outflow rates have been resolved into higher and lower velocity jet 
material. For the clearer case of DG Tau, this revealed that the more collimated 
higher velocity component plays a dominant role in mass and angular momentum 
transport. 
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1. Introduction 



In current star formation theory, jets/outflows from a newly forming star are believed to 
transpo rt significant amounts of energy and momentum away from the region of the central 



source flBallv et all 120071 : iRav et al.l 120071 : IPudritz et all 120071 : IShang et all 120071 ) . This can 



have a big influence on the way in which the stars form, because, for example, jets may 
drive the injection of turbulence in the parent cloud, thereby regulating the star formation 
rate. At the same time, they may be able to extract the excess angular momentum from 
the accretion disk, thus allowing the matter to drift through it and finally accrete onto 
the star. Jets, therefore, are considered a fundamental ingredient in the star formation 
process. To fully understand the mechanisms underlying the physics of these commonly 
observed nebulae, however, it is necessary to know the mass outflow rate, which regulates 
the dynamics of the flow and, therefore, is the most important input parameter for any model 
of flow generation and propagation. For example, in a series of recent papers by our group, 



(Bacciotti et al. 


2002; 


Woitas et al. 


2005; 


Coffev et al. 



2004 ; 2007). The combination of mass outflow rate and toroidal velocity estimates allows 
determination of the angular momentum transported by the flow, and thus also a comparison 
with the angular momentum that the associated disk has to loose in order to accrete at the 
observed rate. Indeed in one test case, namely the jet from RWAur, we found that the 
angular mome ntum present in th e jet is at least 60 - 70 % of that required to be extracted 
from the disk (IWoitas et alll2005l ). This is of course a result of primary importance, since it 
could be the first long-awaited validation of the popular idea that jets exist to remove the 
excess angular momentum in forming system, thus providing a solution to one of the main 
theoretical problems in star formation. 

In order to confirm this finding, however, we have to determine the mass and angular 
momentum outflow rate close to the star, and in a wider sample of objects. A good starting 



1 Based on observations made with the NASA/ESA Hubble Space Telescope, obtained at the Space Tele- 
scope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., 
under NASA contract NAS5-26555. 
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point is to use the HST/STIS spectra that were recently acquired by our group in a survey 
to search for rotation signatures in jets from T Tauri stars. The survey included jets from 
five young stars, namely TH28, DGTau, HH30, CWTau and RWAur. The spectra, taken 
with the slit placed at a few tens of arcseconds from the star, and oriented transversely 
to the flow direction, sampled the collimation region of the jets and allowed us to confirm 
the presen ce of differential D oppler shifts at the jet borders and derive estimates of toroidal 
velocities (jCoffey et al.ll2004l ; 2007). These spectra, taken with the slit perpendicular to the 
jet axis, included several optical forbidden emission lines that can be used to diagnose the 
physical conditions of the gas. Thus, we have embarked in a further analysis of these data 
to understand how the plasma properties behave close to the star where the jet is launched, 
collimated and accelerated, and finally to find the mass and angular momentum outflow rate 
in our sample. We point out that estimating the mass outflow rate at the jet base is more 
likely to exclude spurious effects, such as the addition of mass from possible entrainment of 
ambient gas and dust or significant disruption of the flow, which are less likely to have a 
dominant influence in the early stages of jet propagation. 

The determination of the mass outflow rate requires a knowledge of the total gas den- 
sity, which, contrary to the electron density, is not directly measurable from the observed 
lines. One, therefore, has to use some method to derive the total density from the available 
observational data, and indeed several methods have been developed to this aim in the recent 
past. For example, from a measure of the intrinsic emission line luminosity one can derive 
the number of emitting particles of a given species in the observed volume, which can then 
lead to the gas density under an assumption of abundances. However, this method relies on 
an accurate prior knowledge of reddening estimates, excitation temperature, ionisation state 
of the given species, a nd filling factor, a ll of which bring into the calculation substantial 



uncertainties (see e.g. iNisini et al.ll2005l ). An alternative method is via a determination 



of the hydrogen ionisation fraction of the emitting region. A direct measure of the elec- 
tron density can be easily obtained from the [S II]AA6716,6731 doublet (see e.g. lOsterbrock 
19891 ). Dividing the electron density by the ionisation fraction then closely approximates 
the total gas density. One approach to finding the ionisation fraction is to model the line 
ratios under the assumption of a definite mechanism for the gas heating. There is not yet, 
however, a general consensus on the mechanism causing the jet emission although, without 
a doubt, the observed forbidden lines are excited collisionally. The most widely accepted 
explanation is that the gas is being heated by internal shocks, although other possibilities 
incl ude ambipolar diffusion, turbu lent mixing-layers and compression by jet instabilities (see 
e.g. lLavallev-Fouauet et al.ll2000l and references therein). On the basis that the gas is shock 
heated, lHartigan et al.l (11994 ) constructed a grid of planar shock models to calculate the line 
ratios, and compared the results to spectra of a few stellar jets integrated along the beam. 
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In this way, they were able to find the most likely value of the hydrogen ionisation fraction 
in the flows, and hence the total density, after accounting for the shock compression. This 
method, however, is model dependant, and implies heavy and lenghty calculations. 

Subsequently, a technique was developed whereby the ratios of the optical forbidden 
emission lines of S, N and O in the red wavelength range are used to infer the hydrogen ion- 
isation and the electron temperature without having to assume a specific method of heating, 
but only assuming that the atomic levels are populated co llisionally in the absence of ionising 
photons (IBacciotti et al.lll995t iBacciotti fc Eisloffell Il999l ) . This so-called 'BE technique' is 
based on the recognition that, in the conditions present in stellar jets, the ionisation states 
of oxygen and nitrogen are tightly correlated to that of hydrogen via charge-exchange. Con- 
sequently, the line ratios can be easily modeled and compared with the observed ratios to 
find the most suitable values of the gas physical quantities. In addition to the requirement 
of very little calculations, the technique does not require an accurate knowledge of distance 
or extinction. Furthermore, the results do not depend on any specific mechanism for jet 
formation and/or evolution, and thus can be applied to targets in very different conditions. 



The BE t echnique has been applied to jet observat i ons of varying spectral and spa 
tial re solution (IBacciotti et al.lll995t IBacciotti et al.lll996l ; IBacciotti et al.lll999l ; iNisini et al. 

20061 ). which lately included high angular resolution obs ervations of the 



2005; Podio et al 



DGTau, RWAur, Lk Ha233 jets using HST/STIS multiple-slit configuration, flBacciottill2002 



Melnikov et al.ll2008l ) as well as su b-arcsecond ground-based data of the DG Tau and R W Aur 
jets obta ined with adapt i ve opt ics (lLavalley-Fouquet et al. hood iDougados et al. l boo2h . Very 
recently, lHartigan et al.l (120071 ) used a slightly different version of the BE technique to diag- 
nose the physical condition of the HH30 jet from HST/STIS 'slitless' spectra. In all cases, 
the method has brought very valuable information on the excitation of the jet gas allowing, 
for instance, constraints on models for jet heating, and allowing confirmation that shocks 
are, as expected, the most likely cause of the thermal behaviour of the plasma. In most 
cases, it has also been possible to estimate the mass and momentum flux in the flow from 
the derived total densities. This in turn has facilitated a discourse on the dynamical relation - 
ship between optical jets and coaxial molecular flows (see, e.g. IBacciotti h Eisloffell Il999l ) . 
Furthermore, as mentioned above, in the test case of the RW Aur jet, using the results of the 
BE technique and Doppler gradient measurements, we could give an estimate of the angular 
momentum transported by the flow. 

For all these reasons, we decided to use the BE technique to determine the mass and 
angular momentum flux in our sample of jets observed at high angular resolution. The STIS 
data included all the optical lines necessary to apply this procedure, which has proven to be 
very well suited to the diagnostics of large datasets in extended jet surveys, as is the case 
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of this paper. In the following, we describe the application of the method to our spectra 
(Section[2]), and the obtained results (SectionEJ). Given the form of our input, we present 
the results for the physical conditions as position-velocity (PV) maps, in which the electron 
density, the ionisation fraction, the electron temperature and the total density are resolved 
both as a function of velocity and distance from the jet axis. From these values, we derived 
estimates of the mass and angular momentum outflow rates of the jet close to the launch 
point. Our findings are discussed in Section[U 

2. Observations and Data Analysis 
2.1. Observations 



HST/ST1S spectroscopic observations at optical wavelengths (Proposal ID 9435) were 
made at the base of jets from several T Tauri stars (see Tabled]). These data formed the basis 
of two recently published papers which reported measurements of gradi ents in the Doppler 



shift of the jet radial velocity profile transverse to the flow direction (jCoffey et al.l 120041 ; 
2007). In those papers, the possibility that the measurments can be interpreted as indications 
of jet rotation was investigated. The adopted observing procedure involved centering the 
HST/STIS slit on the T Tauri star, rotating the slit to a position angle perpendicular to 
the jet axis, and then offsetting the slit to a position along the jet which is a fraction of an 
arcsecond from the source, as illustrated in Figured) 

The jets were observed with slit offsets of Of! 3 in the cases of TH28, DGTau, RW Aur, 



but 0'.'2 for CWTau (due to known lack of line emission at 0'/3 (IWoitas et al.ll2002l )). and 
0'/6 for HH30. (The larger offset for HH30 was deemed more appropriate since the star 
and the jet base are obscured by the system's edge-on disk.) These angular distances corre- 
spond to deprojected distances along each jet of 52, 68, 58, 37 and 84 AU respectively. The 
CCD detector was used with the G750M grating, centred on 6581 A, and a slit of aperture 
52x0.1 arcsec 2 . Spectral sampling was 0.554 A pixel -1 , corresponding to a radial velocity 
sampling of ~25 km s -1 pixel -1 , and spatial sampling was 0. "05 pixel -1 . Resolution is twice 
the sampling. 

Long exposures were mad e of individual jet targets, each of ~2000 s to 2700 s (for further 



details see ICoffey et al.l 12004 ; 2007). The data were calibrated through the standard HST 
pipeline, subtraction of any reflected stellar continuum was performed, and signals from 
defective pixels were removed. 

In total, this yielded eight transverse jet spectra, which included Ha, [O I]AA6300,6363, 
[N II]AA6548,6583 and [S II]AA6716,6731. Data analysis is not reported for the faintest 
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Target 


Location 


Distance 

(PC) 


M* 
(Mq) 


Vsys 

(km s —1 ) 


Ijet 
(dog) 


PA iet 
(deg) 


References 


TH 28 


Lupus 3 


170 




+5 


10 


98 


1, 2 


DGTau 


Taurus 


140 


0.67 


+16.5 


52 


226 


3, 4 


CWTau 


Taurus 


140 


1.4 


+14.5 


41 


155 


5, 6, 7 


HH30 


Taurus 


140 


0.45 


+21.5 


1 


33 


8, 9, 10 


RWAur 


Auriga 


140 


1 


+23 


44 


130 


11, 12 



Table 1: Details of TTauri jet targets investigated in this paper. All radial velocity results, 
v rac i, are quoted after correction for the systemic heliocentric radial velocity, v sys . The 
inclination angle of the jet, ij et , is given with respect to the plane of the sky. Values for the 
jet position angle, PAj et , were determined from archival HST images. For our observations, 



Graham & Hever 1988: (2) 


Krautterl '. 


1986; (3) 


Eisloffel et al. 


1998; (4) 


Bacciotti et al. 2002: 


(5) Gomez de Castro 


19931: (6) Hartmann et al. 


1986 


: (7) Hartiean et al. 


200^ 


{ (8) 


Petv et al. 


20061: (9) ADDenzeller et al. 


20051; (10) 


Mundt et al. 


1990; (11 


) Woitas et al. 


2001 


2002; (12) 


Lopez-Martin et al. 


2003. 



Jet Axis 

Receding Jet 



HST/STIS Slit 



TTauri star & disk\ 



Approaching jet 



Fig. 1. — Orientation of the jet and slit. The arrow on the slit indicates the positive direction 
of the y-axis in the position-velocity contour plots. The requested slit position angle was 90° 
with respect to the value of PAj et , Tabled! The actual slit position angles for DGTau and 
HH 30 differ from the requested values by +3° and -6° respectively. This was due to problems 
during observations in finding the right combination of guide stars for the requested angle. 
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targets, namely the approaching jets from TH28, RWAur and LkHa321, due to lack of 
sufficient signal-to-noise for the [N II] emission, which is fundamental to the diagnostic 
procedure used here. 



2.2. Data Analysis 



2.2.1. The BE Technique 
An indepth description of the technique can be found in several p r eviously published 



papers on the meth od and its application (IBacciotti fc Eisloffell Il999l ; iPodio et al.l 12006 
Melnikov et~aP l2008h . Briefly, the BE technique exploits the fact that in low excitation 



conditions, such that the plasma is partially neutral and provided that no strong sources of 
energetic photons are present, the ionisation state of oxygen and nitrogen is dominated by 
charge ex change with hy d rogen . (Justification for this assumption was subsequently con- 
firmed by lHartigan et al.l (120071 ).) In such a plasma, all the sulfur can be considered singly 
ionised. Y oung stellar lets sa tisfy these conditions (if except is made for the so-called irra- 



diated jets iBally et al. 



20071 ). because they have low electron temperatures and effectively 



no photoionisation, as the stellar source is usually a low mass star with little production 
of ionising radiation. Under these assumptions, the ionisation state of both oxygen and 
nitrogen can be expressed as a function of the hydrogen ionisation fraction and the electron 
temperature. Combining this with the fact that the [S II]A6731/A6716 is a well-known in- 
creasing function of the electron density (until it reaches the critical density of [S II], i.e. 
~ 2 10 4 cm -3 ) means that the intensity ratio of any two of these observed emission lines is 
a known function of electron density, hydrogen ionisation fraction and electron temperature. 
In practice, the BE technique consists in calculating a series of model emission line ratios, 
and in comparing the result with the observed ratios, to determine the most suitable values 
of the physical quantities in the emitting gas. In particular, the sulfur doublet ratio is used 
to determine the electron density, while the ionisation fraction and the electron temperature 
are found by combining the ratio of nitrogen to oxygen lines, [N II] /[O I] and the ratio of 
the oxygen to sulfur lines, [O I]/[S II], at the determined electron density ratio. In fact, the 
first ratio is predominantly sensitive to the ionisation fraction ([N II] /[O I] increases with 
the ionisation fraction), while the second is comparably sensitive to both ionisation fraction 
and electron temperature ([O I]/[S II] increases with increasing electron temperature, and 
with decreasing electron density, for densities above the critical density). Once the electron 
density and ionisation are determined, it is possible to calculate the total hydrogen density 
in the emitting gas. 



The application of the technique is straightforward, but it had to be automated for the 



- 8- 



treatment of large datasets. This has been done recently upon the analysis of HST/STIS 
spectra of the LkHa233 jet acquired by our group. A numerical code was written (here- 
after the BE code) which algebraically inverts the functions generated by the line ratios. 
The code was succesf ully tested on the LkHa233 data, and the results were published by 
Melnikov et~all (hoosh . 



We applied this same code to our jet spectra, but we have updated some aspects to 
account for the fact that our data is in a higher excitation regime, i.e. close to the jet 
launch p oint. To com pute the fractional p opula tion in the atomic levels we use a 5-level 
code by iRagal (119921) (see IBacciotti et al.l (119951 ) Appendix A for details). We use solar 
abundances fro m Osterbrock ( 1994T). We update the co llision strengt hs, partially fo l lowin g 
Hartigan et aD (l2007h: O fSlI) fromli ceenan et all Jl996h : fi(N II) from lHudson et al.l (l2005h : 



Q(0 I) from iMendoza et al.l (119831 ). All the other atomic parameters (collisional ionisation 
rates, reco mbination rates, d ielect ronic recombination rates, charge exchange rates) are those 
specified in IBacciotti et al.l (119991 ). 



Before calculating the line ratios, a spurious shi ft in space generat ed by the instrument 
optics was determined for each line and removed (see lCoffey et al.ll2004l ; 2007 for full details). 
Emission lines were resampled in velocity so that differences in dispersion were accounted 
for, and so ratios of emission at the same velocity could be obtained. We considered the 
brightest components in the [O I] and [N II] emission line doublets, the other being always 
one third of the brightest component as dictated by atomic physics. 

In calculating the electron density from the [S II] doublet, values of the [S II]A6731/[S II] A6716 
ratio reaching 2.15 indicate saturation of the doublet ratio such that it is not sensitive 
to higher densities. These positions were thus given a lower limit on electron density of 
2.510 4 cm -3 determined by the [S II] critical density. Meanwhile, values of the same ratio 
barely reaching 0.69 cross the lower limit in density sensitivity, and these positions were 
given an upper limit on electron density of 50 cm -3 . 

In applying the BE code to our jet spectra, we obtained results for electron density, 
electron temperature and ionisation fraction. The electron density was then divided by the 
ionisation fraction to yield an approximation of the total hydrogen density. 



2.2.2. Improvements on the BE code 

While each jet was observed to emit in all forbidden emission lines, it was found that the 
code does not converge at certain pixel positions. The problem was found to arise in cases 
where jet material is travelling with a broad range of velocities, and the various lines peak 
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at different velocity components, i.e. specifically for the DGTau jet. The [S II] emission is 
more intense at lower velocities. Therefore, at higher velocities, the [S II] /[O I] ratio is often 
too low to allow convergence to a single value of the algebraic inversion performed by the 
BE code. 

In an attempt to overcome this problem, therefore, we eased the constraints on the 
algebraic calculations. We take the electron density derived from the [S II] doublet, as before. 
Then, instead of using the BE code, we predict [N II] /[O I] and [O I]/[S II] ratios from a 
mesh of ionisation and temperature values. Comparing the predicted and observed ratios, 
we find the point where the sum of the differences squared is at a minimum. This directs us 
to which combination of temperature and ionisation values best matches the observed line 
ratios. This method of using lighter constraints, which we will call the BE-lite code, produced 
output which was confirmed to be exactly the same as the output from the original BE code, 
when applied to a jet with a narrow range of velocities, i.e. the TH28 jet. 

The BE-lite code was then applied to DG Tau. The original BE code was found to 
represent only 40% of the output of the BE-lite code in the analysis of the DGTau jet. 
The mesh covered a temperature range from 0.210 4 K to 310 4 K in steps of 400 K, and an 
ionisation range from 0.005 to 0.5 in steps of 0.005. 

Upon final application of the appropriate code to each target, the analysis of TH28 
and DGTau gave the best results. This is because these data provide high signal-to-noise 
in all emission lines. The analysis of the other targets was more limited, due to insufficient 
signal-to-noise in [N II] emission, which is less than or equal to that of the other lines in all 
cases. 



3. Results 

All targets examined were found to emit in the forbidden emission line doublets [O I]AA6300,6363, 
[N II]AA6548,6583 and [S II]AA6716,6731. In the adopted observing mode, the slit does not 
collect light from the star. In fact, the slit (of width O'.'l) is placed perpendicular to the 
jet usually at 0'.'3 from the source, and the instrument spatial line spread function has a 
half-width at zero-maximum of 0'.'15 in the optical regime. Therefore, we can be confident 
that the observed emission lines originate in jet material and are not contaminated by stellar 
emission. 

Position-velocity contour maps of the brightest emission line in each doublet are shown 
in Figures (2], [51 [HI [HJ Position- velocity images of the ratios of these lines are shown as 
a guide to the distribution of relative fluxes, Figures El El El E21 CE5J For ease of comparison, 
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ratio maps are presented such that they scale positively with the relavant map of physical pa- 
rameters. (Recall from Section l2.2.1l that [S II]A6731/A6716] increases with electron density, 
[N II] /[O I] increases with ionisation fraction, [O I]/[S II] increases with increasing electron 
temperature and with decreasing electron density) Position- velocity images of jet physical 
parameters are shown in Figures HJ [TOj [13] and [161 These figures represent output from 
the BE code applied to TH28, HH30, CWTau and RWAur, and the BE-lite code applied 
to DG Tau. The typical values are summarised in Table[2j 

The overall results indicate that, at the base of jets of the examined TTauri stars, the 
electron density is high enough to saturate the [S II] doublet (i.e. >210 4 cm -3 ). For all 
targets, the doublet reaches saturation in many points. Where it does not reach saturation, 
the electron density is always at least above 10 4 cm~ 3 . (It should be borne in mind that 
ranges for the other parameters are often derived for regions of saturated [S II] doublet 
ratios.) The typical temperature in this region of the jet is high, and always in the range 
1 - 210 4 K (except for RWAur where it only reaches 510 3 K, although this dataset is the 
poorest). The ionisation level for all targets is in the range 0.03 - 0.3, where DGTau and 
RWAur noticeably define the lower end of this range. The total hydrogen density at the 
jet base lies in the range 5 10 4 -5 10 5 cm -3 , although this is only a lower limit for saturated 
regions. 

Since optically thin forbidden line emission is integrated through the jet, the results 
represent physical conditions seen in projection. This causes a smoothing of any gradients in 
the results. Since such jets are s hown to have an onio n-like structure, i.e. parameters vary 
with each layer from the jet axis faacciotti et aliboooh . the borders of the flow (i.e. regions 



farthest from the jet axis) are less affected by the projection effect, because the line of sight 
does not cut through as many " layers of the onion" . 



3.1. TH 28 receding jet 



The bipolar jet from TH28 was identified by iKrautterl (119861 ). This system is almost 
in the plane of the sky (10°), with a low radial velocity jet (~30kms _1 ), and a narrow 
range of velocities within the jet. Low resolution spectra of this jet have been analysed by 



Bacciotti fc Eisloffell (119991). who find a fall off in electron density with distance along the 



jet already identified by IKrautterl (119861 ). and values for electron density, temperature and 
ionisation in the region of 10 2 cm -3 , 2 10 4 K and 0.07 - 0.61 respectively. These values pertain 
to a distance of 2"-ll" from the driving source. 



Of all our targets, the results for the TH 28 receding jet present the clearest indications of 
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the plasma physical conditions. This is due to good signal-to-noise in all forbidden emission 
lines, comparable jet width in all emission lines, and a similar narrow range of velocities 
in all lines. The position-velocity contour plots, Figured show that the jet emission is not 
symmetric. This is abundantly evident in the [O I] plot, but is also supported by fainter 
[N II] plot. The analysis reveals an asymmetry also in electron density. Although there is 
not much spread in the jet velocity, this asymmetry is clearly evidenced by both lower and 
higher velocity gas. The electron density on one side is a lower limit as it has reached the 
critical density for the [S II] doublet. Aside from the obvious asymmetry, there is no other 
clear trend in electron density, which lies typically in the range of 0.5-1.5 10 4 cm -3 , in the 
region where the sulfur doublet does not indicate saturation. There is, however, a trend in 
the temperature. Values reach 3 10 4 K at a thin front located at the peak red-shifted velocity, 
and then drop to about 2 10 4 K at lower velocities material. The temperature drops further 
still moving from on-axis to the jet borders to 1 10 4 K. Thus, along the jet axis, a spine of 
higher temperatures in the region of 2 10 4 K is decipherable. Also noteable is the decrease in 
temperature with increasing electron density on one side of the jet. This inverse relationship 
is noticeable to a lesser degree over the rest of the jet's cross-section. The ionisation fraction 
reaches 0.5 for higher red-shifted velocities. Otherwise, the ionisation appears constant over 
the whole jet cross-section, with a typical value for this target of 0.3. The combination of 
electron density and ionisation fraction leads to an overall jet hydrogen density of typically 
3 10 4 cm~ 3 . The asymmetry in electron density is reflected giving a lower limit on the 
hydrogen density in a few positions, reaching 1.4 10 5 cm -3 . 

The behaviour of jet physical parameters which we have identified in these high reso- 
lution data s hows some poss i ble ch aracteristics typical of those predicted by resolved shock 



models, e.g. lHartigan et al.l (11994 ) . According to these authors, a resolved shock is repre- 
sented by high temperature which decreases with distance behind the shock, a corresponding 
increase in ionisation, and an increasing electron density as the gas is compressed behind 
the shock. The magnitude of the values depends on the speed of the shock. In our results, 
we see these parameters in velocity space and so if a shock is observed we would expect to 
see a high temperature peak at lower red-shifted velocities. But such a trend is not clearly 
identified, although we so see high ionisation at higher velocities, as would be expected of a 
shock resolved in velocity space. We als o dont see any t r end o f density in velocity space, but 



an asymmetry in the spatial direction. lHartigan et al.l (120071 ) note an agreement with such 



expectations for temperature and ionisation trends in their analysis of low spectral resolu- 
tion data for HH30, but remark on the obvious lack of agreement between the observations 
and models in that the density is expected to increase in the post-shock region but is not 
observed to do so. 



-12- 



3.2. DG Tau approaching jet 



DG Tau was one of the first T Tauri stars to be associated with an optical jet (IMundt &: Fried 



19831 ). This is one of the brightest jets from low mass young stars, and consequently one of 



the best studied amoung young stellar jets. Previous studies of the gas physics, with a view 



to determining the heating mechanism, were conducted by lLavalley-Fouquet et al.l (120001 ) 
based on sub-arcsecond spectro-imaging data taken from the ground. Gas parameters were 
reported as a function of distance from the source and in three velocity intervals. The dat- 
apoints closest to the star are at 0'.'4. Also, Bacciotti et al. (2000; 2002) reported physical 
parameters derived from HST/STIS data of initial portion of this jet with the slit parallel 
to the flow direction. 

The contour plots, Figures[5], reveal that this jet emits in two clearly identifiable velocity 
bins, the higher velocity component (HVC) and the lower velocity component (LVC). [N II] 
is clearly more collimated and traces high velocities, while [S II] is spatially extended and 
traces low velocities. [O I] traces both LVC and HVC with its peak at high velocities, and 
demonstrates that the LVC is indeed less collimated than the HVC. The nominal boundary 
between HVC and LVC was determined to be -lOOkrns -1 , by two-Gaussian deblending of 
the [O I]A6300 line profile. Note that the [O I]A6300 emission is blue-shifted to the edge of 
the CCD. Therefore, for this emission line, we do not detect emission at velocities beyond 
-312kms _1 for instrumental reasons. 

As a direct result of the large velocity dispersion within this jet, and slightly difference 
velocity ranges for different lines, the results for DG Tau presented here are produced with 
the BE-lite code. Immediately, we note the saturation of the [S II] doublet in many points, 
both in the HVC and in the LVC. Because of this saturation, it is not possible to define a 
more detailed trend for electron density as a function of velocity. Spatially, however, we can 
clearly see that th e electron density is h igher close to the jet axis then at the jet borders, 



in agreement with iBacciotti et al.l (120001 ) . A similar spatial trend in other parameters could 
not be established since [N II] was not observed in the jet borders. A typical electron 
temperature of 2.5 10 4 K was found for the HVC and 0.5 10 4 K for the LVC respectively. The 
overall results for ionisation show the observed region of the jet, at Of! 3 from the source, 
to be of low ionisation with values reaching only 0.02 in the HVC and 0.05 for the LVC. 
Recall that we report the electron density in the LVC as a lower limit, implying that the 
temperature and ionisation derived at these positions should be treated with some caution, 
because using the critical density instead of the true density gives high values of temperature 
and low values of ionisation. We then find a total hydrogen density to be quite scattered, 
with a tendancy for the HVC to be more dense. Typical values reach ~ 10 6 cm -3 in the 
HVC and 5 10 5 cm -3 in the LVC. Spatially, no trend in the temperature, ionisation and total 
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density was found. 



Our results agree with trends close to the star reported in the literature. iBacciottil (120021 ) 
report the same low ionisation levels in the same position along the jet as our observations. 
They also find the same trend in the ionisation being higher for the LVC than the HVC, in 
agreement with our results. However, they find higher electron density for the HVC than 
the LVC, whereas we find saturation in both cases. Finally, it is important to highlight that 
the high electron density which w e find at the base of the jet is not found further along the 
jet ( jLavalley-Fouquet et al.l 120001 ) . 



3.3. HH 30 approaching jet 



Mundt fc Friedl (119831 ) identified this jet, upon observing emission line nebulosity in 
the HL Tau complex in Taurus. This jet is in the plane of the sky, with a low radial 
velocity (~ -5 km s -1 ) and n arrow range of velocity within the jet. HST images of this 
system (IBurrows et al.lll996l ) showed for the first time a disk-jet perpendicul ar configuration 



i n you ng s tars. Previous diagno stic studies focused on this target include iBacciotti et al. 



(119991 ) and lHartigan et al.l (120071 ). both of which are based on data of high spatial resolution 
but low spectral resolution. 

There is no division of HVC and LVC in the case of HH30, since the internal velocity 
dispersion is small. Datapoints here are fewer, with respect to the previous two targets, due 
to the low level of [N II] emission. HH30 shows a scatter in the distribution of electron 
density, but with all values above ~10 4 cm~ 3 . Its temperature range is scattered within the 
range 0.4 - 0.810 4 K, and its ionisation lies in the range 0.08 - 0.15. The total hydrogen 
density is thus typically 1.5 10 5 cm" 3 . 



Bacciotti et al.l (119991 ) presented a similar analysis of HH 30, but in one dimension along 
the jet and low spectral resolution, showing how each parameter varies with distance from 
the source. Close to the source, they report an electron temperature of 2 10 4 K, and a lower 
limit on electron density due to saturation of the [S II] doublet ratio. They report that the 
jet ionisation fraction rapidly rises from 0.065 at 0'/2 to 0.1 at Of! A, and then slowly increases 
up to 0.140 within 2" from the source. Our results are in good agreement with these values 
for the examined position. 



More recently, lHartigan et al.l (120071 ) conducted a study of the HH30 jet using HST 
'slitless' spectroscopy. This method relies on the fact that choosing a wide (2") slit ensures 
the jet is unresolved spectrally, but produces sufficient dispersion to separate images of the 
jet in the different optical emission line. The data represents the jet to 4"(600 AU) along the 
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flow, and in two epochs. Hence, the analysis shows how jet parameters vary with distance 
and time. Our results are in excellent agreement with those reported at 0'.'3 from the source 
in both epochs (i.e. log electron density ~4.6, electron temperature ~5000K, ionisation 
fraction ~0.1). 



3.4. CW Tau approaching jet 



The jet from this source was discovered in optical images by iGomez de Castrol 11993 
Although this j et is well studied both through high resolution imaging and spe ctroscopic ob- 
servations (e.g. iDougados et al.ll2000l ; lHartigan et al.ll2004l ; ICoffey et al.l 120071 ). no previous 
analysis of the gas physics has been published for this target. 

CWTau very clearly exhibits a separation of jet emission into a HVC and LVC, with 
the nominal boundary between the two measured as -50kms _1 using two-Gaussian deblend- 
ing. As with HH30, points are few due to the level of [N II] emission. From the available 
datapoints, both velocity components appear saturated in electron density with some unsat- 
urated points in the LVC at typically 0.8 10 4 cm -3 . Values are found in electron temperature 
of 0.810 4 K and ionisation was typically of 0.2 - 0.3. Total hydrogen density then became 
810 4 cm" 3 . 



3.5. RW Aur receding jet 



RW Aur is a complex triple star system. The highly collimated jet from RW Aur A was 
i dentified by iHirth et al.l (I1994J ) . Previous studies of the gas physics have been conducted 
( IDougados et al.l |2002| ). and gas parameters reported as a function of distance from the 
source and in three velocity intervals. In that study, the datapoints closest to the star are 
at 0"4. 

The analysis of the RW Aur jet yields just a few datapoints in the diagnostic maps, due 
to the poor [N II] emission. The [S II] doublet begins to saturate and leaves a scattered 
pattern which can be vaguely observed to increase with increasing red-shifted velocities. 
The derived temperature reaches 0.5 10 4 K. Like the DGTau jet, RWAur also indicates low 
ionisation at the observed location a long the jet, with a m aximum of 0.07. Our results are 



in agreement with previous studies (IDougados et al.l 120021 ) . which give upper limits at the 
jet base to be ~ 10 4 K for the temperature, ~ 10 -2 for the ionisation fraction, ~10 4 cm -3 for 
the electron density. The resulting total hydrogen density reaches 5 10 5 cm -3 . 
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3.6. Errors 



Sources of error on the estimated parameters depend on the flux measurement error, 
the effects due to differential extinction, the choice of elemental abundances. 

The error on all fluxes is calculated as the rms of the background, after cleaning cosmic 
rays and defective pixels. This varies according to target, but remains relatively constant 
across all emission lines for a given target spectrum. From this we obtain the three sigma 
threshold reported as the contour floor in each contour plot of the emission lines. Errors 
on each of the derived physical parameters due to uncertainty on the line flux measurement 
were found to be on the order of 30% for points of good signal-to-noise. 

The amount of dust extinction towards the various positions along these jets is unknown , 
and may even vary along the beam of a single object (see lNisini et al.ll2005l ; iPodio et al.ll2006l ). 
The extinction towards the targets in this study could not be calculated, since the emission 
lines observed are too close to each other in wavelength. Therefore, we have not accounted 
for the effects of extinction by dereddening our emission line fluxes. Differential reddening, 
i.e. the variation of extinction with wavelength, does affect flux ratios. However, since we are 
dealing with emission lines spanning a small range in wavelength, the effect is not expected 



to be significant. Previous studies confirm that this is indeed the case. lHartigan et al. 



( 120071 ) examined the effects of differential reddening on their study of HH30. They find 
that the uncertainty introduced by an extinction of A„~ 1 is not significant compared to 
the uncer tainties in the relative abu ndances, and hence choose to ignore the effects in their 
analysis. iBacciotti &: Eisloffell (119991 ) report that, assuming a fiducial extinction of Aj,~3, 
variations in relative flux values reveal errors in the final jet parameters of amounts no larger 
than the measurement errors, being at most 10% for the ionisation frac tion and 15% for 
the t e mperature. For each of our sources, the level of extinction is A,,<3 (jGraham &: Hever 
19881 : iBeckwith et~aD Il99pl : lHartigan et~al~1l2007l : lEiroa et alihopj bhez et al.lll997h . Fur- 
thermore, the visual extinction of T Tauri stars is expected to sharply decrease moving away 
from the source, and hence should be much lower than 3 at 0'.'3 - 0'.'6. Therefore, we neglect 
the error introduced by the fact that we are not dereddening the fluxes prior to the analysis. 

Lastly, the BE code relies on an assumption of elemental abundances for the determina- 
tion of temperature and ionisation. Solar abundances are ad opted, but may not accurately 
reflect conditions in star forming clouds. IPodio et al.l (120061 ) show that, using abundances 
determined for the interstellar medium in Orion, the values of ionisation and temperature 
inferred using the BE technique are within 15% of those obtained assuming the most recent 
determinations of solar abundances. Since our errors arising from flux measurements are 
in the region of 30%, we determine that only an insignificant error is introduced to our re- 
sults by adopting solar abundances, instead of the actual abundances for the considered star 
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forming regions (i.e. Taurus-Aurigae and Lupus). 



4. Discussion 

As discussed in SectionHJ the fundamental parameter to be derived from the diagnostic 
analysis is the mass outflow rate of the flow. Mass flux was determined from an average of the 
total hydrogen density, n(H), the poloidal velocity of the emission peak for [O I], v p , and the 
jet radius, rj et , in order to get an average value integrated over the full jet cross-section. The 
formula used was Mj et = fim p n(H)iTr'j et Vp, where /i is the mean molecular weight (~ 1.24) 
and m p is the mass of a proton. The jet poloidal velocity was determined from radial 
velocity measurements as v p = v ra d/ sin i j Pt . For HH30, this value is taken from proper 



motion measurements close to the source (IBurrows et al.lll996l ). since the inclination angle 



is so small that poloidal velocities cannot be determined from radial velocity measurements 
without incorporating considerable error. The jet cross-section was determined by assuming 
that the FWHM of the [OI]A6300 emission can be taken as the jet diameter, since it traces 
both HVC and LVC In cases where the jet material is clearly divided into HVC and LVC, 
i.e. for DGTau and CWTau, the mass outflow rate was estimated for each component 
separately. The estimates are reported in Tabled Since the measurement uncertainties on 
the quantities in the calculation are at the 30% level, the uncertainty on the calculation of 
the mass outflow rate leads to a factor of a few, similar to previously published estimations. 
We note that any mass outflow rate estimate conducted in this way can represent an upper 
limit since, not only does this method assumes a homogenous distribution of material within 
the jet, but also the densest regions dominate the observed emission leading to a bias towards 
higher densities. However, close to the star we do expect a homogenity and high density. 
Conversely, the positions where the [S II] doublet saturates indicate a lower limit on the 



Target 


n e 


T e 


X £ 


n(H) 




(10 4 cm" 3 ) 


(10 4 K) 




(10 4 cm" 3 ) 


TH 28 receding jet 


1.0 


2.0 


0.30 


3 


DG Tau approaching jet HVC 


>2.5 


2.5 


0.02 


90 


DG Tau approaching jet LVC 


>2.5 


0.5 


0.05 


50 


HH30 approaching jet 


1.5 


0.6 


0.1 


15 


CW Tau approaching jet HVC 


>2.5 


0.8 


0.30 


8 


CW Tau approaching jet LVC 


0.8 


1.2 


0.20 


4 


RW Aur receding jet 


>2.5 


0.5 


0.06 


30 



Table 2: Summary of typical values of the output from the BE technique analysis, as shown 
in detail in FiguresHJ [TD1 [13] and (THl The lower limit on electron density is not carried 
through to the total density for reasons explained in SectionHJ 
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Fig. 2. — Position-velocity diagrams for the brightest component in each of the three 
forbidden line doublets, for the TH28 receding jet. Contours begin at the three-sigma 
level of 1.35 10~ 15 and are in steps of 10% up to the peak levels of 6.36, 3.55 and 
2.01 10 -14 erg cm -2 s _1 A -1 arcsec 2 for each plot respectively, while the two-sigma level is 
marked by the dotted line. 
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Fig. 3. — Position- velocity flux ratio images demonstrating the useable points in each ratio, 
for the TH28 receding jet. The plots are calculated based on only the brightest line in the 
[O I] and [N II] doublets, from which the second component of the doublet was inferred using 
atomic physics. Positions of zero value indicate pixels with emission flux below three sigma 
for either one or both species in the ratio. 
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Fig. 4. — Position- velocity images showing physical jet parameters as a function of distance 
from the jet axis and velocity, for the TH28 receding jet. The plots represent electron 
density, electron temperature, ionisation fraction and hydrogen density calculated for each 
pixel of the data. 
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Fig. 5. — Same as Figure[2]but for the DGTau approaching jet. Contour floor is 2.7 10~ 15 
and ceilings are 42.8, 7.9 and 15.5 10~ 14 ergcm -2 s _1 A^ 1 arcsec 2 . Note that the [O I]A6300 
emission is blue-shifted to the edge of the CCD. Therefore, for this emission line, we do not 
detect emission at velocities beyond -312 km s -1 for instrumental reasons. 
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Fig. 6. — Same as Figure[3]but for the DGTau approaching jet. 
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Fig. 7. — Same as FigureHJbut for the DGTau approaching jet. 
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Fig. 8. — Same as Figured but for the HH30 approaching jet. Contour floor is 2.4 10 
and ceilings are 10.9, 1.7 and 9.2 10 -14 erg cm -2 s _1 A^ 1 arcsec 2 . 
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Fig. 9. — Same as Figure[3]but for the HH30 approaching jet. 
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Fig. 10. — Same as FigureHJbut for the HH30 approaching jet. 
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Fig. 11. — Same as Figure[2]but for the CWTau approaching jet. Contour floor is 2.10 10 15 
and ceilings are 2.8, 1.2 and 1.2 10~ 14 ergcm -2 s _1 A^ 1 arcsec 2 . 
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Fig. 12. — Same as Figured but for the CWTau approaching jet. 
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Fig. 13. — Same as FigurefTOl but for the CWTau approaching jet. 
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Fig. 14. — Same as Figured but for the RW Aur receding jet. Contour floor is 3.0 10 15 and 
ceilings are 16.3, 0.8 and 1.8 10~ 14 erg cm -2 s _1 A -1 arcsec 2 . 
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Fig. 15. — Same as Figure[3] but for the RW Aur receding jet. 



-26- 



0.6- 
u 0.4- 

P 

a 

« 0.2- 
< 

3 0.0- 

E 

o 

8 -0.2- 

C 

PJ 

Q -0.4- 
-0.6- 



RW Aur Receding Jet 



Electron Density (10 4 cm ~ 3 ) 



Electron Temperature (10 4 K) 



0.4 

0.3 

0.2 

0.1 
0.0 



h 



-200 



200 



-200 



200 



0.6 



0.4 



0.2 



0.0 



E 
o 

8 -0.2 

C 

B 

S -0.4 
-0.6 

Fig. 16. 



Ionisation Fraction 



0.06 

0.05 

0.04 

0.03 

0.02 

0.01 
0.00 



n 



Hydrogen Density (10 4 cm ~ 3 ) 



-200 200 -200 200 

Radial Velocity (km s" ) Radial Velocity (km s" 1 ) 

Same as FigurefTOl but for the RWAur receding jet. 



40 




30 




20 




10 










-27- 



electron density and hence on the total density measurement. Therefore, considering these 
counter-acting factors, the mass outflow rate is more of an average indication, rather than 
an upper or lower limit. 

Angular momentum flux estimates were obtained by assuming the radial velocities and 
mass outflow rates were indicative of gas physics over the full jet cross-section. Values 
were determined by combining the mass outflow rate with toroidal velocities derived from 
the differences in radial velocity across the jet. The a verage toroidal ve locity of the jet is 
calculated from the transverse radial velocity gradient (jCoffey et al.ll2004l ; 2007) and the jet 
inclination angle as = (Av ra d/2)/ cos ij e t- The radius from the jet axis was taken as half 
the FWHM, as measured for the mass outflow rate calculation. Since we approximate the 
mass outflow rate and toroidal velocity, and presume them to be unvarying over the full cross- 
sectional area of the jet, we are satisfied with a simplified formula for the integration of the 
angular momentum outflow rate, Lj P ± = -\I i, it',,' 4 i, t- F or HH30, no radial velocity difference 
were reported above the error bars (jCoffey et al.ll2007l ). and so no angular momentum outflow 
rate was calculated. The results are reported in Table[3j 



4.1. TH 28 receding jet 

For this jet, the mass and angular momentum outflow rates in this lobe are of the 
order of those previously reported for TTauri jets, although on the low side, with values of 
1.2 lO~ 8 M yr _1 and 2.5 10~ 6 M Q yr^ 1 AUkms -1 respectively. The asymmetry in electron 
density identified in this target translates to a mass outflow rate which is doubled on one side 
of the jet axis compared with the other side. Values were calculated for two semi-circular 
areas (seen in projection) either side of the central pixel row which represents the jet axis. For 
the calculation, the poloidal velocity of the emission peak was used, along with the average 
hydrogen density either side of the axis. This yielded estimates of 0.4 and 0.9 10~ 8 M Q yr _1 
for the two semi-circles. Furthermore, the latter value relies heavily on regions where the 
electron density is high enough to cause saturation of the [S II] doublet ratio and so it is 
likely to represent a lower limit (if the jet is homogenous at this position). 



4.2. DG Tau approaching jet 

For DG Tau, the jet material is clearly divided into HVC and LVC, and so the mass out- 
flow rate was estimated for each component separately. The mass outflow r ate is higher in 



the HVC at 4.1 lO" 8 M yr- 1 compared to 2.6 lO^M^r" 1 in the LVC. iBacciotti et al. 
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( 120021 ) take a nominal boundary of -70kms _1 between the medium velocity component 
and LVC (instead of -100 km s" 1 ), and see the trend reversed. They find values of 5.1 
and 8.3 lO~ 8 M yr _1 as determined from HST/STIS spectra with the slit placed paral- 
lel to the flow axis. Our total mass outflow rate of 6.7 10~ 8 M Q yr _1 is somewhat lower 
than 1.3 10 _7 M Q yr _1 for the parallel slit data, but the latter assumed jet widths for 
each velocity component which were based on slit positions, and which transpired to be 
1.5 times our F WHM measurements. The mass accretion rate for this system is measured as 
2 10 _6 M o yr _1 (IHartigan et al.lll995l ) and so, doubling our values to account for a bipolar 
flow, we find an mass flow ratio of Mj et /M acc ~ 0.07. Our results show that the jet angu- 
lar momentum outflow rate in the approaching jet is also considerably higher in the HVC. 
We found estimates of 8.4 and 4.5 10 -6 M yr _1 AU kms -1 for HVC and LVC respectively. 
The total of 1.3 10~ 5 M Q yr _1 AU kms -1 represents the highest angular momentum outflow 
rate of all our target s. The value is in good agreement with 3.8 10~ 5 M yr _1 AUkms -1 
(IBacciotti et al.ll2002l ). considering the aforementioned difference in the adopted jet radii. 
Our high spatial resolution results show higher values of angular momentum outflow rate 
closer to the axis, which is traced by the HVC. In this jet HVC and LVC are well separated 
and so projection effects may be lessened, given that the HVC traces the core of the jet while 
the LVC traces the outer borders. 



4.3. HH 30 approaching jet 



The mass outflow rate of HH30 is very low at 4.0 10 -9 M yr _1 , and it lies outside the 
typical range for T Tauri stars of 10~ 6 to 10 -8 M yr _1 . Nevertheless, our result compares well 
with a previous d etermination from Plateau de Bure interferometric data of 1.0 10~ 9 M yr _1 
(jPety et al.l 120061 ) . Our results are also in ag reement with the esti mate of 1.7 10 -9 M yr _1 , 
made assuming a jet velocity of 200 km s -1 (IBacciotti et al.lll999l ). No angular momentum 
outflow rate could be calculated for this target due to lack of evidence of toroidal velocities 
JCoffev et al.ll2007h . 



4.4. CW Tau approaching jet 



As in the case of the DGTau jet, HVC and LVC are clearly separated in CWTau jet. 
The mass outflow rate is found to be almost evenly distributed between the two components 
giving values of 0.4 and 0.3 lO _8 M yr _1 for HVC and LVC respectively but here, unlike 
DG Tau, the datapoints are few and so the result is less rel iable. When compared with 
the mass accretion rate of lO~ 6 M yr _1 (IHartigan et al.lll995l ). we obtain a ratio of ~0.01. 



-29- 



As in all cases where the surfer doublet ratio saturates, the mass outflow rate of the jet 
may be underestimated. The angular momentum outflow rate in the approaching jet, also 
evenly distributed in velocity, was found to be 0.5 and 0.6 10 -6 M Q yr _1 AUkms -1 for the 
HVC and LVC respectively. The total of 1.1 10~ 6 M Q yr _1 AUkms -1 represents the lowest 
angular momentum outflow rate of all our targets. 



4.5. RW Aur receding jet 



Despite the lack of datapoints, a mass outflow rate of 1.7 10 -8 M Q yr -1 was estimated 
at 0'.'3 above the disk-plane. This value is lower than 3.0 10 _8 M Q yr _1 reported at ft' 2 
above the disk-plan e from HST/STIS spectra with the slit parallel to the flow direction 
( jWoitas et al.l 120021 ). As in the case of DGTau, the difference arises from the fact that 
the latter calculation adopted a jet radius according to the slit position rather than the jet 
FWHM. Their jet radius is 1.5 times our measurement. Our mass outflow rate result, when 
doubled to account fo r a bipolar flow, is 2% of the mass accretion flux of 1.6 10~ 6 M Q yr -1 
( IHartigan et al.lll995l ). We estimate the angular momentum outflow rate in the receding jet 
to be 2.9 10 -6 M yr -1 AU km s _1 , compared to 1.0 10 ~ 5 M yr -1 AU km s _1 estimated for the 
aforementioned parallel slit data (jWoitas et al.ll2002l ). Again, the difference arises from the 
difference in radii. The estimates increase to 6.0 10~ 8 yr" 1 and 2.6~ 5 M ( nvr~ 1 AUkms -1 
for measurements further above the disk-plane at 0'.'5 (jWoitas et al.ll2005l ). 



Target 


FWHM 

(arcsec) 


v p 

(km s _1 ) 


Mjet 
(10- 8 Mgyr" 1 ) 


v$ 
(kms -1 ) 


Ljet 

(10" 6 M yr" 1 AU kms" 1 ) 


TH 28 receding jet 


0.32 


165 


1.2 


8 


2.5 


DG Tau approaching jet HVC 


0.14 


224 


1.1 


20 


8.4 


DG Tau approaching jet LVC 


0.20 


97 


2.6 


12 


4.5 


HH30 approaching jet 


0.20 


54 


0.4 






CW Tau approaching jet HVC 


0.17 


148 


0.4 


10 


0.5 


CW Tau approaching jet LVC 


0.35 


42 


0.3 


7 


0.6 


RWAur receding jet 


0.18 


129 


1.7 


14 


2.9 



Table 3: Mass, Mj et , and angular momentum outflow rates, Lj et , for each jet target. Mass 
flux is integrated over the full jet cross-section. Values are divided into higher velocity 
component (HVC), and lower velocity component (LVC) where applicable. 
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5. Conclusions 

We have analysed the gas physics for several TTauri jets close to the launching point, 
based on high resolution HST/STIS spectra taken with the slit perpendicular to the flow 
direction. We have applied the BE diagnostic technique to the line emission spectra, to 
obtain maps of the electron density, temperature and ionisation fraction. To this aim, we 
have also adapted the BE code, to allow inclusion of more datapoints in the calculations. The 
new approach involves easing the algebraic constraints, which cause difficulties where jets 
exhibit a broad range of velocities and the various the various lines trace different velocity 
components, i.e. DGTau. Our results represent the first survey of physical conditions at 
the base of T Tauri jets presented in the form of position-velocity diagrams for the physical 
quantities. In fact, with our high spatial and spectral resolution dataset, we resolved the jet 
physics as a function of velocity and distance from the jet axis, in the region just a few tens 
of AU above the disk plane where the flow is launched. 

The overall survey results indicate that, at the jet base, the plasma has a high electron 
density (>2 10 4 cm -3 ), high electron temperature (210 4 K), and low ionisation level (0.03 - 
0.3) which varies considerably depending on the target. Indeed, in all cases saturation of 
the [S II] doublet is reached in at least some of the datapoints, and so in this region close 
to the launch point of the jet the electron density reported is often a lower limit. This is 



different to findin gs further along the jet which are a factor of ten lower (IBacciotti et al.lll999 



Podio et al.ll2006l ). In the case where previous studies have been carried out, we find good 
agreement for values reported close to the star. In the case of TH 28 (which presents the best 
dataset), possible shock signatures are present, thus providing a observational indications 
that shocks can contribute to heating the jet close to the source. 

We determine the mass and angular momentum outflow rates for the jets close to their 
base. Estimates we determined for the mass and angular momentum outflow rates, both 
of which are fundamental parameters in constraining models of accretion/ejection struc- 
tures, particularily if the parameters can be determined close to the jet footpoint. Val- 
ues for a single jet lobe are in the range 4.0 10~ 9 to 6.7 10 _8 M Q yr _1 and 1.1 10~ 6 to 
1.3 10~ 5 M Q yr _1 AUkms -1 . Again we find good agreement with the literature in cases where 
values were reported close to the star. Mass flow ratios were found to be Mj et /M acc ~ 0.01 - 
0.07 where accretion rates were available in the lit erature (i.e. for 3 of 5 targets). This is in 



the range predicted b y accretion-ejection models (IKonigl &: Pudrita |2000| ; ICasse &: Ferriera 



200dlShu et al.ll2000h . 



Although we have examined a region of the jet at about 50 - 80 AU from the source 
corresponding to the collimation zone, we note that the region where the jet is formed and 
launched is believed to be on scales of less than 1 AU, a region which is currently out of 
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the reach of present instrumentation and often obscured by infalling matter. We await 
near-infrared interferometry as an opportunity to observing this zone. 
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